clear;clc;

load v_dz_2.mat;

Q=9.73*(v/1000).^(2.85);

[zz,xx]=size(v);

dx=12.5;
dz=4.0;
x_tick=[0:100:xx];x_label=dx*x_tick;
z_tick=[0:100:zz];z_label=dz*z_tick;
Font=14;

figure;
set(gcf,'units','normalized','position',[0.1 0.1 1.0 0.5]);
subplot(1,2,1);
imagesc(v);
ylim([0,zz]);xlim([0,xx]);
set(gca,'xaxislocation','top');
set(gca,'tickdir','out');
xlabel('x(m)','FontSize',Font);
ylabel('Depth(m)','FontSize',Font);
set(gca,'FontSize',Font);
set(gca, 'XTick', x_tick);
set(gca,'XTickLabel',x_label,'FontSize',Font);
set(gca, 'YTick', z_tick,'FontSize',Font);
set(gca,'YTickLabel',z_label);

subplot(1,2,2);
imagesc(Q);
ylim([0,zz]);xlim([0,xx]);
set(gca,'xaxislocation','top');
set(gca,'tickdir','out');
xlabel('x(m)','FontSize',Font);
ylabel('Depth(m)','FontSize',Font);
set(gca,'FontSize',Font);
set(gca, 'XTick', x_tick);
set(gca,'XTickLabel',x_label,'FontSize',Font);
set(gca, 'YTick', z_tick,'FontSize',Font);
set(gca,'YTickLabel',z_label);